Local Superfluidity at the Nanoscale 
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We have performed quantum Monte Carlo simulations measuring the finite size and temperature 
superfluid response of 4 He to the linear and rotational motion of the walls of a nanopore. Within the 
two-fluid model, the portion of the normal liquid dragged along with the boundaries is dependent 
on the type of motion and the resulting anisotropic superfluid density saturates far below unity 
at T = 0.5 K. The origin of the saturation is uncovered by computing the spatial distribution of 
superfluidity, with only the core of the nanopore exhibiting any evidence of phase coherence. The 
superfluid core displays scaling behavior consistent with Luttinger liquid theory, thereby providing 
an experimental test for the emergence of a one dimensional quantum liquid. 
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Superfluidity or dissipation-free flow, is rooted in quan- 
tum mechanics with the wave function of the entire fluid 
being described by an emergent global macroscopic phase 
0. In 4 He, this breaking of gauge symmetry has dra- 
matic consequences for the liquid below the superfluid 
transition, T\ ~ 2.17 K. It is well established that su- 
perfluid helium can flow through extremely narrow con- 
strictions, impenetrable to the normal liquid, with a ve- 
locity v s = (h/m) J V6 limited only by a critical velocity 
first understood by Landau. In "rotating bucket" ex- 
periments, where a container of superfluid is rotated at 
an angular frequency lj, vortices can be spontaneously 
created, yielding a non-zero quantum of fluid circulation 
k — § v s ■ dr = —W, where W G Z is the topolog- 
ical winding number, equal to the number of vortices 
within a closed circulation loop. The quantitative details 
of the superfluid state were first probed in the celebrated 
Andronikashvili torsional oscillator experiment in 1946 
where it was determined that a superfluid fraction of the 
total fluid does not contribute to the classical moment 
of inertia. This observation led to the development of 
Tisza's phenomenological two-fluid model where the su- 
perfluid state is understood as two intertwined liquids, 
having normal (p n ) and superfluid (p s ) components with 
total density p = p n + p s . 

The superfluid-normal transition in bulk 4 He is in the 
three dimensional (3d) XY universality class. As the spa- 
tial dimension of the system is reduced, the enhance- 
ment of fluctuations suppresses the transition tempera- 
ture to zero in 2d and precludes any long range ordered 
state in one dimension. It is intriguing to consider how 
the continuum two-fluid picture holds up in the low di- 
mensional limit where the bosonic helium system should 
be described by the universal harmonic Luttinger liquid 
(LL) theory pQ at low energies. Such a correlated liquid 
is strongly fluctuating with any phase coherence decay- 
ing algebraically as a function of distance at T = K. 
Experimental realizations of low dimensional bosonic sys- 
tems have been achieved in ultra-cold atomic gases [2] at 



low densities where the interactions are expected to be 
weak and short ranged, but at higher densities, direct ob- 
servations of LL behavior are still lacking. Previous ex- 
perimental work in low-<i quantum fluids has focused on 
superfluid helium confined to porous materials with a ra- 
dial length scale in the nanometer range, and there is evi- 
dence for a new quantum phase occurring at low tempera- 
ture |3, 4 . Recently, experiments have demonstrated the 
feasibility of measuring the superflow of helium through 
nanometer sized holes [5 . The next step will be to sys- 
tematically decrease the radius of the nanopore, thereby 
providing a quasi- Id flow geometry in which the super- 
fluid properties of helium can be measured. Unequivocal 
evidence of LL behavior in 4 He filled nanopores will re- 
quire a detailed understanding of the signatures of low 
dimensional superfluidity in the crossover regime and we 
have performed large scale numerical simulations to elu- 
cidate this behavior. After providing details of our mi- 
croscopic model and numerical methods, we present the 
low temperature superfluid density which we identify as 
originating from only the Luttinger liquid core region of 
pores with nanometer radii. 

The starting point is a system of helium-4 confined in- 
side a nanopore of radius R and length L formed as a 
cylindrical cavity in a slab of amorphous silicon nitride. 
The interactions between helium atoms are modeled via 
the Aziz potential [6] while confinement is achieved by 
combining the effects of short range repulsion with the 
walls of the pore and a long range dispersion force be- 
tween helium and the surrounding medium [7 j. The re- 
sulting quantum many-body Hamiltonian can be exactly 
simulated using continuous space Worm Algorithm (WA) 
quantum Monte Carlo (QMC) [8]. Within the path in- 
tegral formulation, this method exploits the quantum- 
classical isomorphism, performing Metropolis sampling 
of (d + l)-dimensional configurations of bosons that can 
be visualized as worldlines or trajectories in an imagi- 
nary time (r) direction [9]. For helium at finite temper- 
ature, the worldlines obey a periodicity condition in the 
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FIG. 1. (color online). The origin of superfluidity in 4 He 
filled nanopores measured by the winding (W, left column) 
and area (A, right column) estimators in the path integral 
representation. At high temperature, helium worldlines are 
short, containing only single particles and there is no super- 
fluid fraction: p s /p = 0. As the temperature is lowered, the 
worldlines may link with each other, winding around periodic 
boundary conditions in the simulation cell and having a large 
projected net area perpendicular to the axis of rotation. 



additional dimension modulo identical particle permuta- 
tions, owing to their bosonic symmetry. The superfluid 
response of the low temperature system is directly linked 
to the existence and properties of long connected world- 
line exchange cycles consisting of many individual atoms. 

Our simulations employ a fixed chemical potential 
/i//uB = —7.2 K to ensure helium atoms in pores of 
L = 75 A and R = 3.0 — 15.0 A are in thermal con- 
tact with a bath held at saturated vapor pressure (SVP) 
for temperatures between 0.5 — 2.25 K. 

In QMC, the superfluid density is measured using lin- 
ear response theory by considering the effects of bound- 
ary motion [10]. Within the two-fluid model, it is sup- 
posed that a superfluid fraction p s /p will remain station- 
ary while the normal portion p n / ' p will be dragged along 
with the walls of the container. In the nanopore geome- 
try, two types of motion, depicted in the first row of Fig. [I] 
are possible: translational motion of the walls along the 
cylindrical axis and a rotation around it. These two types 
of response correspond to different measurements in the 
QMC, related to the geometry and topology of particle 
worldlines, represented as closed loops due to periodicity 
in imaginary time. In the case of transverse wall motion, 
(left column, Fig. [TJ we consider a container with peri- 
odic boundary conditions along the axis of the pore, and 
for the purposes of visualization, imagine the volume in- 
side the pore to be mapped onto the surface of a torus. 



The major circumference of the torus is equal to L while 
the minor one is hj3 with f3 = l/k^T. In the path inte- 
gral representation, W = (l/£)X^i Jo^ dr[dzi(r)/dr] 
where Zi(r) is the z-component of r^(r), the (d + 1)- 
dimensional position of particle z, is equivalent to the 
number of times the imaginary time trajectories of the N 
particles wrap around the periodic boundary conditions 
of the sample. In WA simulations, the imaginary time r 
must be discretized and we use kB^r/h = 0.004 K _1 to 
minimize Trotter error. The resulting superfluid density 
pY is related to the variance of the distribution of wind- 
ing numbers present [9] [10] through a winding estimator, 
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where m is the mass of a helium atom and (• • • ) indi- 
cates a QMC average. At high temperature, the helium 
atoms behave classically, with spatially localized world- 
lines containing only a single atom. Long exchange cycles 
are extremely unlikely and so (W 2 ) = 0. As the temper- 
ature is lowered, the kinetic energy can be reduced by 
linking worldlines together. Such particle exchanges can 
be efficiently sampled within the WA using spatially local 
updates, producing configurations with extended world- 
lines that wind around the periodic boundary conditions 
((W 2 ) 0), producing a finite superfluid response. 

An alternative approach, equivalent in the d > 3 
thermodynamic limit [TT], measures the non-classical re- 
sponse of the fluid to a small rotation. The superfluid 
fraction p^ / p is then equal to the non-classical rotational 
moment of inertia fraction (I c \ — I)/ I c \ where / is the ob- 
served moment of inertia and I c \ is the total classical 
moment of inertia. The superfluid fraction defined in 
this way can be estimated in the QMC by measuring the 
worldline area A of closed particle trajectories projected 
onto a plane perpendicular to the axis of rotation [T2] 
through an area estimator, 
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For a rotation about the z-axis: A = 
amElJo^dHriir) x [dr^/dr]},. At high 
temperature, the projected mean squared areas are 
uncorrelated and y/ (A 2 ) ~ A 2 where A = y/2nh 2 j3 / m 
is the thermal de Broglie wavelength and the resulting 
superfluid fraction will be vanishingly small for large 
pores as seen in the right column of Fig. [I] As the 
temperature is reduced and long exchange cycles become 
energetically favorable, there is a distribution of finite 
projected areas and p^/p > 0. The caveat to this ap- 
proach is that angular momentum conservation requires 
that the dimensions of the simulation cell perpendicular 
to the axis of rotation do not have periodic boundary 
conditions. A quantum fluid confined inside a nanopore 
thus provides a geometry where we can directly compare 
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FIG. 2. (color online). A comparison of the superfluid fraction 
of helium confined inside nanopores with L — 75 A (symbols) 
measured using the winding number and projected area of 
particle worldlines with the experimentally measured value 
[13] (line) for bulk 4 He at saturated vapor pressure. Upper 
right insets show instantaneous quantum Monte Carlo particle 
configurations projected into the z — plane while the lower 
left inset in the first cell details superfluidity on a finer scale. 
All panels share the legend shown in the first row. 




FIG. 3. (color online). The confined superfluid core seen 
by comparing the particle p, winding superfluid and area 
superfluid number densities measured via quantum Monte 
Carlo simulations for a R = 13 A radius nanopore of length 
L — 75 A at T — 0.75 K. The first row shows the axially 
averaged local densities inside the pore while the symbol plot 
contain an additional angular average. 



the two estimators of the superfluid fraction as shown 
in Fig. [2] The symbols correspond to QMC measure- 
ments performed using Eqs. and Q and the insets 
in the upper right corners are instantaneous particle 
configurations projected on the plane z = 0. The solid 
line is the experimentally measured superfluid fraction 
of 4 He at SVP taken from Brooks and Donnelly [13] for 
comparison. For pores with R < 9 A we do not observe 
any persistent superfluid response above T = 0.5 K 
while for R > 10 A, pY ,A /p becomes non-zero at a 
radius dependent onset temperature shifted below the 
bulk value of T\. This is the expected behavior for 
a quantum fluid constrained inside a porous material 
where confinement has been shown to reduce T c and 

Ps/P- 

The observed superfluid response measured via the 
winding number is effectively one dimensional, originat- 
ing from flow along the pore axis, and any non-zero value 
should be considered a finite size effect that will disap- 
pear as L — )> 00. Another striking feature distinguishing 
the nanopore superfluid fraction from that of bulk he- 
lium is an apparent saturation at low temperature for 
both pf 1 and at a value much less than one. A hint at 
the origin of this behavior can be observed by examining 
the spatial configurations inside the pore in Fig. [2] The 
interplay of interactions between helium atoms as well 
as with the surrounding amorphous SiaN4 leads to states 
exhibiting a series of cylindrical shells [14] |15] , equivalent 
to the formation of thin film layers of helium observed on 



2d substrates including silicate [T6] . 

The competition between the tendency of bosons to 
delocalize at low T and the strong geometrical confine- 
ment effects in the pore can be investigated by measuring 
the local contribution of the superfluid density. This was 
achieved by histogramming the radial r-dependence of 
the winding number [17] or path area [18] 

nY(r)="^<W-W{r)) (3) 

with W and A the full pore values defined above 
while W{r) = {2TrrL 2 )- 1 Y.f =l .Q f3 dT[dzi{T)/dT}5{r - 
ri{r)) and A(r) = (AirrL)- 1 £ti ^ dr^r) x 
[dri(r)/dr]} z S(r — r^{r)) where I c \(r) — mr 2 is the 
classical moment inertia of a single helium atom and 
rj-{r) = ^x 2 (t) + y 2 {r). A comparison of the average 
particle number density with the two types of local su- 
perfluid density can be seen in Fig. [3] for a nanopore with 
R = 13 A, L = 75 A at T = 0.75 K. The upper pan- 
els have been averaged over the z-axis while the lower 
plot contains the fully cylindrically symmetric radial val- 
ues defined in Eqs. ([3| and Q. The saturation of the 
total superfluid density seen in Fig [2] can now be imme- 
diately understood in terms of a spatial "phase" sepa- 
ration where only the inner volume of the pore contains 
superfluid helium while the outer shells remain classical, 
adhering to the walls. The results are qualitatively sim- 
ilar for R = 10 — 15 A with the two outermost shells 
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making a negligible contribution to the superfluid den- 
sity. The local superfluid estimators are nearly iden- 
tical while their total values can differ by a factor of 
two. This presents no paradox due to their different 
normalizations: p^ / ' p = (2ttL/N) rdrp™ (r), p^/p = 

(2irLm/I c \) rdr[r 2 pf(r)] needed to account for local 
contributions to the classical moment of inertia present 
in inhomogeneous fluids I c \ = 27rLra rdr[r 2 p(r)]. 

The core region exhibiting a non-zero superfluid re- 
sponse is nearly one dimensional, having a radius of 
R < 6 A. In d = 1, fluctuations preclude the exis- 
tence of any long range superfluid order, and instead, 
the helium system should be described at lowest or- 
der, by the linear quantum hydrodynamics of Luttinger 
liquid (LL) theory pQ with effective Hamiltonian H = 
{hv/2ir) J L dz[(d z <l)) 2 /K+K{d z 0) 2 ]. The phases </>(z) and 
0(z) are defined in terms of the second quantized he- 
lium field operator ^(z) ~ ^d z 6(z)e l ^ z ^ such that 
[4>(z),d z '9(z f )] = i7rS(z — z'). Its low energy modes have 
dispersion e(k) — hvk and the value of the Luttinger pa- 
rameter K tunes the system between algebraic superfluid 
(K <C 1) or solid (K ^> 1) order. For a real physical 
system, the velocities vj = v/K and vn = vK can be 
related to the microscopic parameters of the underling 
many-body Hamiltonian. By comparing the predictions 
of harmonic LL theory, derived from the grand partition 
function Z = Tr exp[— f3(H — pN)) with the measure- 
ments from finite temperature QMC simulations, vj and 
Vn can be determined. For quasi- Id helium confined in- 
side nanopores with R < 3 A, this has already been ac- 
complished [19] but for larger radius pores, required the 
use of an ad hoc cutoff radius when analyzing QMC data. 
The physical origin of this cutoff is now fully understood 
as the radius of the superfluid core and we expect it to 
be described by LL theory [20 : 
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where Os(z^q) is a Jacobi theta function with f ^(z^q) = 
d 2 0s(z,q) and p c = (N c /N)p where 7V C is the number of 
atoms in the core. For each radius, we have performed 
a rescaling of the total superfluid response displayed in 
Fig. [2] and determined the velocity vj(R) through fitting, 
that yields the best collapse of all low temperature data 
onto Eq. ([5| . The results are displayed in Fig. [4] where 
the temperature scaling of the nanopore superfluidity is 
consistent with Luttinger liquid theory. 

Much remains to be done, including a confirmation of 
the predicted pore length scaling of pf 1 / p. In addition, 
it seems natural to contemplate the effects of disorder, 
surely present in the pore walls, as well as the introduc- 
tion of fermionic 3 He which may strongly alter superflu- 
idity as bosonic exchanges will be suppressed in Id. 

In conclusion, for pores with R = 1 — 1.5 nm we 
have observed a finite and anisotropic superfluid response 
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FIG. 4. (color online). The superfluid fraction of the core of 
nanopores for varying radii which can be collapsed onto the 
universal prediction from Luttinger liquid theory. The inset 
shows the extracted value of the phase velocity hvj obtained 
by fitting to the winding number estimator for each radius. 



above T — 0.5 K, with a magnitude that is dependent 
on whether linear or rotational motion of the cylinder is 
considered. The difference is large and arises from the 
absence of any classical moment of inertia in the truly 
Id limit where flow is still possible. Experiments prob- 
ing this breakdown of the naive two-fluid picture could be 
performed by comparing the superfluid fraction measured 
by capillary flow and a nanoscale Andronikashvili tor- 
sional oscillator. Our results also indicate that when the 
radii of nanopores becomes sufficiently small, the super- 
fluid fraction may exhibit plateaus, increasing in steps, 
due to the classical sticking of wetting layers near the 
pore walls. This is in stark contrast to the usual smooth 
temperature dependence of p s /p observed for bulk 4 He 
and could provide a signature of the crossover to Id be- 
havior. If the fraction of atoms adhering to the nanopores 
walls could be discerned, possibly by comparing flow 
rates at high and low temperature, an examination of 
the finite size and temperature scaling of the superfluid 
density could confirm that confined low-dimensional he- 
lium is a Luttinger liquid. 
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